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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 
RESEARCH MEMORANDUM 


A THEORETICAL ANALYSIS OF THE DISTORTION 
OF FUEL-SPRAY-PARTICLE PATHS IN A 
HELICOPTER RAM-JET ENGINE DUE 

TO CENTRIFUGAL EFFECTS 


By S. Katzoff and Samuel L. Smith, III 
SUMMARY 


Certain anomalies in the behavior of a ram-jet engine tip-mounted 
on a rotating helicopter blade had led to the suggestion that centrifugal 
effects were distorting the fuel-spray pattern. The present paper is a 
theoretical study of these effects. The differential equations of motion 
of the spray particles were set up and integrated numerically for sixty- 
four selected combinations of the various parameters. For greater ease 
of computation, all but one of these cases were confined to two-dimensional 
motion of the particles (that is, motion in the plane of the rotor). A 
simple method for determining the approximate final direction of travel 
of the particles, without solving the equations of motion, was developed. 


Plots of the solutions show that centrifugal action does have a 
marked effect in deflecting the larger particles toward the outer wall 
(that is, the wall of the combustion chamber that is farthest from the 
center of the rotor) and thus out of the main combustion region. Reduc- 
tion of the particle size dlone probably would not solve the problen, 
since the fine particles may not move far enough laterally from the 
nozzle before being carried downstream in a rich-mixture cylindrical 
core near the center of the combustion chamber. Several suggestions 
for alleviating the fuel-distribution problem are given. 


INTRODUCTION 


In some recent experimental studies of a ram-jet engine mounted on 
the tip of a helicopter blade, it appeared that the thrust obtained was 
less than that measured in bench tests of the same engine in an air stream 


RESTRICTED 


2 RESTRICTED NACA RM L53A02 


of velocity equal to the rotor tip speed; also, the combustion was more 
erratic. One suggested explanation of the difference is that the fuel 
particles are thrown radially outward by centrifugal action, with 
resulting asymmetric fuel distribution and accumulation of fuel toward 
the outer wall (that is, the wall of the combustion chamber that is 
farthest from the center of the rotor). Contour-like stains on the 
outer wall of the jet unit, presumably caused by a pool of liquid fuel, 
served to substantiate this theory. 


The present analysis was undertaken to determine the plausibility 
of this theory, to determine the range of fuel particle sizes for which 


such an effect would be appreciable, and to learn the relative importance 


of some of the other parameters involved, Essentially, the analysis 
consists in setting up the differential equations of motion of the fuel 
droplets in a rotating flow field and integrating for several values of 
some of the major parameters. The air flow in the ram jet was assumed 
to be nonturbulent and moving &long circular streamlines with uniform 
angular velocity about the center of the rotor. Because of these and 
other idealizations, the results, which are presented as plots of fuel- 
particle paths, cannot be considered as quantitative representations of 
actual paths; the results do, however, define the range of parameters 
within which centrifugal effects become important. : a 


For further simplification of the computations, the majority of 
the cases considered were confined to two-dimensional motion of the fuel 
particles; that is, motion in a plane perpendicular to the axis of rota- 
tion. Sixty-three two-dimensional cases were calculated, corresponding 
to selected combinations of three particle sizes, two rotor radii, two 
angular velocities, two internal-flow velocities, two fuel injection 
velocities, and three fuel injection angles (all in the upstream direc- 
tion). One three-dimensional case, corresponding to injection of the 
fuel particle out of the plane of rotation of the nozzle, was calculated 
to give a general indication of the effect of a third (normal) component 
of the initial velocity on the particle path. n. 


SYMBOLS = 
X, y, Z rectangular coordinates in fixed reference frame; origin at 


center of rotor (see fig. 1) 


5,05 2 cylindrical coordinates, relative to reference frame rotating 
&bout z-axis at center of the rotor (see fig. 1) 


r, Y, 2 coordinates of the fuel particles in a rotating reference 
frame moving with the fuel nozzle (see fig. 1) 
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w angular velocity of rotation of the r, 0, z system with 
respect to the x, y, z system 

(o) angular velocity of rotation of the spray nozzle relative to 
the r, 0, z system 

Pa air density 

Pe fuel density li 

v kinematic viscosity of air 

D diameter of spherical fuel particle 

R Reynolds number of fuel particle 

CD drag coefficient of fuel particle 

B initial injection angle of fuel particle, measured with respect 
to the forward direction (normal to radius), positive outboard 

t time 

m mass of the fuel particle 

g acceleration of gravity 

Subscripts: 


O 


indicates initial value (at the time of ejection) 
ANALYSIS 


tions.- The analysis was made on the basis of the following 


simplifying assumptions: 


(1) The fuel particles are spherical, so that available data for 


the drag coefficient of a sphere may be used 


(2) Evaporation of the droplets is negligible; that is, the radius 


of the droplet is constant 


(3) The effect of gravitational force on the particles is negligible 


(4) The air flowing over the fuel spray is without turbulence and 


has a constant velocity with respect to the nozzle 
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Ihe first three assumptions are considered to be satisfactorily valid 

over the region of travel of the fuel particles prior to burning. The 

fourth assumption is much more arbitrary but is considered satisfactory 

for the present purpose. ;  /. JJ.  »;PwmD 4.0 m 
Sphere drag.- It is necessary to determine the drag of a sphere as 

a function of size and velocity in order to set up and solve the equa- 

tions of motion of the fuel particles. A considerable simplification 

of the equations of motion would have resulted had it been permissible 

to use Stokes' classical law of drag of a sphere (because of the simple 

proportionality between drag and velocity). Unfortunately, however, 

the Reynolds numbers encountered (20 <R < 1000) extend well beyond the 

range for which Stokes' law or any of its several modifications is valid, 

and empirical data, such as that given in references 1, 2, and 3, had to 

be used. The empirical curve of drag coefficient against Reynolds number 

from reference 1 is shown as a dotted line for this range in figure 2. 

Since the equations of motion were to be solved numerically on the Bell 

Telephone Laboratories X-66744 relay computer at the Langley Laboratory, 

an algebraic expression for the drag was required. The empirical drag 

curve was approximated in the desired interval by the equation 


72.56 


D = Rx 13.0k 13.0% + 0.415 (1) 


for purposes of computation. The curve of this equation is also shown 
in figure 2. Its maximum error in the range 20 < R < 1000 is about - 
5 percent and its mean error is less than 3 percent. Actually the 
initial Reynolds numbers for a few cases were somewhat above 1000, but 
the resulting error was small, that is, well within the accuracy of the 
analysis. 


The velocity that produces drag is, of course, the particle veloc- 
ity relative to the local (internal) air. In the following analysis 
the three mutually perpendicular components of this relative velocity 
are TP, ró, and z. Therefore, the resultant relative velocity will 


be Jr? + r262 + 22, which, together with the particle size and the 
kinematic viscosity of the air, determines the Reynolds number and 
hence the drag coefficient. The net drag, which is directed so as to 
oppose the relative motion of the particle, is given by 


P 2 e > e. 
Drag - = x Cp (+2 + rege + 22) E (2) 
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and the three mutually perpendicular components of the drag are given 
by 


2 acl. . š : 
(Drag), = an Cp Ee + 126 y 32 (2a) 
(Drag) o = au Op rÓjr? + rèdè + 2 (2b) 
Pa «Dé fe ; r 
(Drag), = 2 Cp 2 FÈ + 7262 + 22 (2c) 


Differential equations of motion of spray particles.- For con- 
venience the motion of the spray particles will be considered relative 


to a reference frame rotating with an angular velocity equal to the 
angular velocity of the internal air in the ram-jet unit. For example, 
for a ram-jet unit operating at a tangential velocity (tip speed) of 
800 fps with respect to air at rest and an internal air-flow velocity 
of 200 fps relative to the nozzle and walls, the internal air has an 
absolute velocity (with respect to a fixed reference frame) of 600 fps, 
so that the angular velocity of the rotating reference frame has the 
value w mou radians per second, where ro is the radial distance 
O 
to the center of the fuel nozzle; therefore, in the general case, 


_ Tip speed - Internal-flow velocity 
Initial radius : 

reference frame presents some computational advantage, since the veloc- 

ity of the particle relative to this system is the velocity that produces 

drag. Furthermore, since these drag forces tend to constrain the par- 

ticle to move with this flow, it is this value of w that determines 

the eventual centrifugal acceleration of the particle. 


W The choice of this rotating 


In this rotating reference frame, let the particle coordinates be 
r, €, and z so that, relative to a stationary reference frame 
(x, y, z), the coordinates are r, O0 +t, and z, where t is time 
elapsed since these two frames coincided. The following relations, 
therefore, connect the rotating cylindrical and stationary rectangular 
systems: 


r cos(8 + wt) 


Pa 
Il 


y =r sin (6 + wt) 


Zz =Z 
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Differentiating twice to obtain the acceleration gives: 


F cos(0 + wt) - 24(6 + m)sin(e + wt) - rid +w) cos(0 + wt) - 


xX = 
ró sin(0 + wt) 

y = f sin(o + wt) + 2r(6 + w)cos(o + wt) - r( + o)?sin(O0 + wt) + 
rö cos(0 + wt) 

7 sg - 


In order to find the radial and tangential components of the 
acceleration at time t, assume that at that instant the radius coincides 
with the x-axis, so that X and y will be the two desired components. 
At this instant (0 + wt) is zero and the equations reduce to 


Sa . 2 
Y - r(0 +0) 


X = 
y = rò + 2r(Ó +0) 
Z =Z 


Thus, the radial acceleration consists of Y plus the centrifugal 
acceleration, and the tangential acceleration consists of r0 plus the 
Coriolis acceleration. These three acceleration components are now 
related (since force = mass x acceleration) to the previously derived 
drag components (eqs. (2a), (2b), and (2c)) to give the following 


equations: 
Pa (io EYES + r262 + 22 -— r(ò + o) 
"zu /0 P Spes 
2 [1 e a vw a » 
- pa (9 róyr + r282 + ZÈ = YO + 2r(8 + w) 
Pa, apé AE 252 "a T T 
- rx uw Cp zyr^ + red tZ = Z 
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where 


72.56 


Cp "Ua 1 y + 0.415 


and where the mass of the particle m is expressed as 


B xD? pe 
6 


m 


Substituting for m and rearranging gives: 


Po a Ch iyi? + r202 + 22 - r(Ó + wé = 0 (3a) 


rö + Pa Cp róyr? + r262 + z^ + ar(6 + w) =0 (3b) 
kp pd 
x Cy GE? + 120% + 22 = 0 (3e) 


These three equations constitute the system of differential equations 
that must be solved in order to obtain the paths of the particles. For 
the cases of two-dimensional motion of the particle (z constant), equa- 
tion (3c) is eliminated and the term 22 is dropped from the remaining 
two equations. 


Z + 


Solution of equations.- As was mentioned previously, most of the 
calculations were performed on the Bell relay computer. The process 
was essentlally & step-by-step numerical integration of the simultaneous 
equations (3a) and (3b) by Heun's method (ref. 4). For the calculation 
of any particular case, w is given, as previously described, by the 
rotor tip speed, the internal airspeed relative to the walls of the jet 
unit, and the radius ro from the rotor axis to the fuel nozzle. The 
assumed outward and upward components of the fuel-injection velocity 
relative to the nozzle give the initial values of f and 2. If the 
forward injection-velocity component relative to the nozzle is given by, 
Say, vi, then the total tangential velocity of the particle relative 
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to the fixed frame is the tip speed + vi, which must be the initial 
value of (W + O)ro. That is, initially 


Tip speed + vt = (w + do)ro 


Thus 
due (Tip speed - wro) + vi 
O To 

or 


Internal flow speed (relative to wall) + vy 
r 


O 


The initial values of 0 and z may be taken as zero. For each path 
the selected initial conditions were specified and a suitable interval 
for the integration steps was chosen. In all the two-dimensional cases, 
time increments of 0.0004 second were used for the early parts of the 
paths, and increments of 0.0008 second were used after the particle had 
approximately established its final direction. The one three-dimensional 
calculation was performed manually by Euler's method (ref. 4), with 
0.0002-second time increments for the early parts and 0.0004-second 
increments for the later parts of the path. 


These calculations give the particle paths in the original rotating 
coordinate system (r, @, z). Before plotting the particle paths, however, 
it was considered desirable to transform to a system moving with the ram 
jet itself, so that the curves would show the paths as they would appear 
to an observer riding on the ram-jet nozzle. This transformation may be 
easily performed by shifting the angular displacement of each point of 
the path by ot, where o is the angular velocity of the ram jet rela- 
tive to that of the internal air in the ram jet, and t is the time 
required for the particle to travel to that point from the injection 
nozzle. This new angular location is that given in figures 3 to 19 
and is designated by Y, where Y is related to € by the equation 


y=0 - ot 
RESULTS 
Cases calculated.- The values of the parameters and initial condi- 


tions used to calculate the sixty-three two-dimensional cases considered 
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in the present report are presented in table I. The first thirty-six 
cases listed in this table are all possible combinations of three 

particle sizes, three injection angles, two rotor radii, two tip speeds, 
one injection velocity, and one internal-flow speed. The remaining 
twenty-seven cases were considered to ascertain the effect of increasing 
the internal-flow velocity or the fuel injection velocity or both. The 
centrifugal accelerations at the fuel nozzle for the cases calculated 

are 13308 for ro = 15 feet and 800 fps tip speed, 795g for ro = 25 feet 


and 800 fps tip speed, 2070g for rg = 15 feet and 1000 fps tip speed, 
and 1240g for rg = 25 feet and 1000 fps tip speed. The cases covered 


are considered to be fairly representative of the region of present 
interest and to define the limits of the problem. 


The fuel density was constant at 1.61 slugs/cu ft throughout all the 
calculations. The values of the air constants used were: air density, 


0.00238 slug/cu ft; and air kinematic viscosity, 0.000172 ft?/sec. These 
values of the air characteristics correspond to operation at altitudes 

of about 4000 feet for the 800 fps tip speed, and 9000 feet for the 

1000 fps tip speed. These altitudes are only approximate since they 
depend on inlet losses and atmospheric conditions. 


The values for the case of three-dimensional motion are: particle 
size D = 0.00% inch; rotor radius ro = 25 feet; tip speed, 800 fps; 
internal-air velocity with respect to the nozzle, 200 fps; fuel-injection 
velocity relative to the nozzle, 150 fps; injection direction, in the 
vertical plane through the ram-jet axis normal to the radius (p = 0°) 
but inclined 45% from the plane of the rotor. 


Results of the calculations.- Figures 3 to 18 present the particle 
paths for the sixty-three two-dimensional cases that were calculated, 
and figure 19 presents the particle path for the three-dimensional case. 
It should be noted that the curve corresponding to f = 0°, D = 0.00} inch, 
is given on both figures 7 and 8 for ease of comparison. This repetition 
also occurs on figures 11 and 12 and on figures 15 and 16. ‘The outline 
of the walls and flame holders of the particular ram-jet unit that occa- 
sioned this investigation is dotted in on figure 3(a). It will be noted 
that the plotted paths extend well into the combustion chamber. Actually, 
these later portions of the paths can not be considered generally valid, 
because of combustion, turbulence and other effects; also, many particles 
may be caught on the flame holders. These portions are shown here, how- 
ever, not only for the sake of completeness, but also because they might 
be applicable for larger ram-jet units, with greater distance between 
the nozzles and the flame holders. 


The calculated points sre shown on the paths in order to indicate 


the time history of the particles and to provide means for estimating 
velocities and accelerations. For the two-dimensional cases the points 
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represent time intervals of 0.0004 second. over the first, highly curved, 
parts of the path, and 0.0008 second over. the later, nearly straight 
parts. The point where the interval is doubled is easily identifiable 
on each path. Similarly, the projection lines in the three-dimensional 


plot of figure 19 indicate time intervals of 0.0002 second for the initial | 


portion of the curve and 0.0004 second for the later portions. 


DISCUSSION 


Particle-path distortion.- All the curves (figs. 3 to 19) show 
similar characteristics. There is sn initial sharply curved part, in 
which the general direction of travel changes from forward to rearward, 
followed by the approach to a fairly well-defined final direction. The 
distortion of the particle paths due to centrifugal force is clearly 
evident on all the curves. The centrifugal forces effectively super- 


impose an outward motion to the particles, so that the curves are deviated 


from the center and have an outwardly inclined final direction. Increased 
tip speed and decrdased rotor radius, both of which result in increased 


centrifugal effects, are accompanied by increased displacement and inclina- 


tion of the path. In this connection, a comparison of the curves of fig- 
ure 3 with the corresponding curves of figure 6 reveals that the paths 

are nearly identical, especialiy for the smaller particles; a simple 
calculation shows that the centrifugal acceleration (ar) is almost the 
same for both cases, even though the radii and tip speeds are quite dif- 
ferent. Increased internal-flow velocity (compare, for example, fig. 11 
with fig. h(b)) results in less centrifugal distortion since the particles 
are turned and carried downstream more rapidly; thus, the centrifugal 
forces have less time in which to act. Increasing particle size has the 
opposite effect (see, for example, fig. 8) - since the mass (and hence 


MOET) increases much more rapidly with particle diameter than does the 
air drag, the centrifugal distortion is accentuated. — ~ 


With the lower internal air speed, only the paths of the smallest 
particles undergo deviations sufficiently small to appear reasonably 
acceptable. These smallest particles, however, do not penetrate very 
far outward before they are turned downstream (see fig. 3(a) or 5(a)), 
so that use of a very fine spray may result merely in & narrow column 
of very rich mixture passing through the combustion chamber. Increased 
ejection velocity of the fine particles somewhat increases the width of 
the column and thereby helps alleviate the problem; however, the use of 
several fine-spray nozzles distributed across the duct area ahead of the 
combustion chamber'is probably the most direct solution. Increasing 
both the particle size and the internal air speed (figs. 15 to 18) also 
results in a broader column without excessive outward inclination of the 
paths; these higher internal air speeds id: eee a probably too high 
for stable burning. 
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Although the present results clearly indicate the nature of the 
problem imposed by the centrifugal effects, it is apparent that the 
analysis was considerably idealized; and some judgment, based on an 
understanding of the nozzle spray characteristics (particle size and 
velocity distribution) and of the flow within the ram jet, would be 
required in applying the present results in the development of a heli- 
copter ram jet. For example, the internal-flow velocity is not uniform, 
as has been here assumed, but is high in the inlet just ahead of the 
nozzles, decreases as it passes down the diffuser to the combustion 
chamber, and then increases again as combustion proceeds along the 
chamber. In addition there is a low-velocity stagnation region imme- 
diately ahead of the nozzle and its attachment brackets, so that the 
outward penetration of the particles immediately after ejection may be 
much more than here shown; this effect would result in increasing the 
width of the fuel spray colum that passes into the combustion chamber. 
Furthermore, instead of all particles being ejected from one point, as 
assumed in the present calculations, there is frequently a cluster of 
nozzles, the spread of which adds to the width of the fuel spray column. 
Finally, it should be pointed out that the implied assumption of a smooth 
laminar type of flow past the nozzles and through the combustion chamber 
is considerably at variance with the actual state of the flow. The 
large-scale turbulence produced by the flame holders will serve to broaden 
the distribution of the spray particles. The distribution of the finest 
particles will be broadened considerably beyond that indicated in, say, 
figure 3(a) and the turbulence may even suffice to prevent some or most 
of the intermediate particles from being thrown directly to the outer 
wall by the centrifugal effects. It seems doubtful, however, that the 
turbulence could similarly keep an appreciable proportion of the larger 
particles from being thrown to the wall. In this connection an addi- 
tional effect should be noted. Because the centrifugal acceleration of 
the air within the ram jet may be of the order of 500 to 1000g, the 
radial pressure gradient within the ram jet may correspondingly be 
500 to 1000 times the vertical pressure gradient of the atmosphere. 
Rough calculations show that, in such pressure gradients, the air-density 
variations resulting from burning can result in convection currents 
comparable in violence to the turbulence introduced by the flame holders. 


With regard to general methods of alleviating the problem presented 
by centrifugal effects, the use of fine-spray nozzles, high ejection 
speeds, and distributed arrangment of the nozzles have already been 
indicated.  Ejecting most of the fuel inward and grouping the nozzles 
inboard of the ram-jet axis should also be helpful. Another possibility 
is skewing the inlet and combustion chamber so that the walls would be 
more nearly parallel to the skewed (or outwardly inclined) particle 
paths. Theoretically, there need be no corresponding loss in thrust for 
such a design as long as the exit direction remains normal to the rotor 
blade; however, the external drag of such a configuration may be exces- 
Sive. In any case it seems likely that an adequate solution may be 
obtained for any given design problem by following these or similar 
approaches. 
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Limiting direction of particle paths.- It is interesting to note 


that the later portions of the three fuel particle paths corresponding 
to the three injection angles appear to he approximately parallel; that 
is, their directions are approximately independent of the initial injec- 
tion angle. Since this direction is fairly indicative of whether the 
particle will pass into the combustion zone or be lost on the outer 
wall, an effort was made to develop a simple method for approximating 
this direction without actually going through a solution of the dif- 
ferential equations. Such a method would..be useful for a rough evalua- 
tion of the path for cases that are not within the range covered e 
figures 3 to 18 of the present report. - 


In & strict sense, & "final" direction may not be clearly definable - 
the particle path (relative to the nozzle) probably approached some 
spiral asymptote, the exact shape of which is somewhat irregular owing, 
for example, to the variation of Reynolds number along the path. In any 
case, the present problem is not concerned with the final direction _ 
when r approaches infinity, but rather with the direction when r is 
of the order of rg (or perhaps a few inches more than ro, as indicated 
on the curves of figures 3 to 18) and the path curvature has decreased 
practically to zero. Specifically, it is proposed to put the second 
derivatives equal to zero in the differential equations (eqs. (3a) 
and (3b)), put r equal to rg (or perhaps a few inches more than ro), 
and solve the two equations algebraically for r and 6, which, together 
with r, determine the direction of the particle. 


Letting Y =6 =0 in equations 3(a) and 3(b) gives 


ac 
iz 


Division of equation (la) by (hb) gives S. ji 


od Hê + r208 - aa + OP = ) | (ka) 


5 2 + reò? + 2 + 1262 + 2(è + IT = | | (hb) 
272 = .r?0(0 + w) (5) 
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Solving equation (5) for r, and substituting into equation (lb) gives, 
after rearrangement and simplification 


(ae - W) = -2(8 + w) (-alè + 0) 
p 


Squaring both sides and rearranging gives 


; 9 
è — h(9 t c) (6) 


Equation (6) can be solved directly for the desired value of a; however, 
the fact that Cp is also a function of Ò would make the computation 
quite laborious. Since the desired final value of 6 is small rela- 
tive to w, the equation can be simplified for computational purposes, 
without undue loss of accuracy, to 


kwé 
2,2 3Pa, Er 
D LoD 


Me 
tl 


or 


MEN 
3Pa 
vi) 


where the need for choosing the negative square root is apparent from 
equation (5). It must be remembered that Cp appearing in equation (7) 


is not a constant, but a function of the velocity (through the Reynolds 
number) and thus of @. 


Me 


(7) 


The Reynolds number may be found from the equation 


= Dr [6 (è - W) 
ars unm 
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is the total realtive velocity in terms of 0, (derived 


by substituting from equation (5) into the expression Vr? + r8? ) and 

may be found by referring to either figure 2 or equation (1), pro- . 
vided that R is in the valid range, 20 <R < 1000, which appears to be 
the range of interest for the present problem. Equation (7) may be 
readily solved for the value of ð by a sintple iteration method (ref. 5). 
From the value of # thus derived, the value of f -may then be obtained _ 
by use of equation (5); thus 


where r AE SE 


6(6 + w) i 
-2 : -(8) 


These values of è and rt together with the radius will determine the 
final direction in the frame moving with the internal air. The final 
direction of the particle travel in the reference freme Sona with the 
ram-jet unit is given by the angle Ed x = 


or 


where the latter form does not require the intermediate calculation of r. 
This method was checked for several cases against the paths calculated by 
the Bell computer. The inaccuracy was found to be small (about 3 percent) 
for the smallest particles (D = 0.002 in.) but to increase to about _ 

15 percent for the largest particles (D = 0.006 in.). Thus, the method 
is probably satisfactorily accurate for the range of basic parameters, 
particularly particle size, of practical interest. 


An approximate relation between the final radial velocity r and : 
the centrifugal acceleration wèr may be derived by substituting for 8 


from equation (7) into equation (8). Equation (8) then becomes yt des s 
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from which, since 6 is small relative to w, it is apparent that 


r x ufr 


That is, the final radial velocity is approximately proportional to the 
square root of the centrifugal acceleration. 


CONCLUDING REMARKS 


Simultaneous differential equations, amenable to numerical integra- 
tion, have been developed for the fuel-spray-particle paths in a ram-jet 
engine operating on the tip of a helicopter rotor blade. A method of 
determining the approximate asymptotic direction of the motion has also 
been given; that is, a method of finding the approximate outward devia- 
tion of the path in the combustion chamber. 


Solutions of the equations for a range of conditions show that the 
centrifugal effect tends to throw the heavier fuel particles to the 
outer wall of the ram-jet unit, while particles fine enough to minimize 
this effect may not reach far enough laterally from the nozzle to provide 
a favorable fuel distribution in the combustion chamber. Several methods 
of overcoming the problem are suggested, such as using a distribution of 
fine-spray nozzles and skewing the combustion chamber outward. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va. 
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TABLE 1.- TABLE OF PARAMETERS AND INITIAL VALUES FOR CASES CALCULATED 
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SCalculations were made for all possible combinations of the values in this line 
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Figure 1.- Sketch of the relationship between the three coordinate systems. 
The z-axis is normal to plane of paper. 
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Figure 2.- Experimental sphere-drag coefficient against Reynolds numbers, 
and comparison curve for algebraic expression used in calculations. 
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Figure 3.- Paths of fuel particles ejected in plane 
0%, end 45°. Rotor radius (to nozzle), 15 feet; 


internal air speed rel 


ative to nozzle, 200 fps; 


relative to nozzle, 150 fps; centrifugal acceler 


of rotor, with B = -150, 
tip speed, 800 fps; 
fuel-injection velocity 
ation at nozzle, 1330g. 
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(a) Particle diameter, D, 0.002 inch. 


Figure }.- Paths of fuel particle 


s ejected in plane of rotor, with È = -A59, 


09, and 45°. Rotor radius (to nozzle), 15 feet; tip speed, 1000 fps; 
internal air speed relative to nozzle, 200 fps; fuel-injection velocity 
relative to nozzle, 150 fps; centrifugal acceleration at nozzle, 20708. 
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Figure l.- Continued. 
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Figure 4.- Concluded. 
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(a) Particle diameter, D, 0.002 inch. 


Figure 5.- Paths of fuel particles ejected in plane of rotor, with f = 459, 
0°, and 45°. Rotor radius (to nozzle), 25 feet; tip speed, 800 fps; 
internal air speed relative to nozzle, 200 fps; fuel-injection velocity 
relative to nozzle, 150 fps; centrifugal acceleration at nozzle, 795g. 
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Figure 5.- Continued. 
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(c) Particle diameter, D, 0.006 inch. 


Figure 5.- Concluded. 
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(a) Particle diameter, D, 0.002 inch. 


Figure 6.- Paths of fuel particles ejected’ in plane of rotor, with B = -459, 
0%, and 45°. Rotor radius (to nozzle), 25 feet; tip speed, 1000 fpa; 
internal air speed relative to nozzle, 200 fps; fuel~injection velocity 
relative to nozzle, 150 fps; centrifugal acceleration at nozzle, 1240p. 
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Figure 6.- Contimed, 
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Figure 7.- Paths of 0.004-inch-d 
of rotor, with 8 = -h59, 09, and 45°. 


Radial distanos from rotor axis, r, 


diameter .fuel 


ed NARA 16. 


particles ejected tn pfane of 
Rotor radius (to nozzle), 


15 feet; tip speed, 1000 fps; internal air speed relative to nozzle, 
200 fps; fuel-injection velocity relative to nozzle, 250 fps; centrif. 


ugal acceleration at nozzle, 20708. 
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Figure 8.- Paths of 0.002-, 0.00%-, and 0.006-inch diameter fuel particles 


ejected in plane of rotor, with B = 09. Rotor radius (to nozzle), 

15 feet; tip speed, 1000 fps; internal air speed relative to nozzle, 
200 fps; fuel-injection velocity relative to nozzle, 250 fps; centrif- 
ugal acceleration at nozzle, 2070g. 
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Figure 9.- Paths of 0.004-inch-diameter fuel particles ejected in plane 
of rotor, with B = -45° and 45°. Rotor radius (to nozzle), 25 feet; 
tip speed, 800 fps; internal air speed relative to nozzle, 200 fps; 
fuel-injection velocity relative to nozzle, 250 fps; centrifugal 
acceleration at nozzle, 795g. 
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Figure 10.- Paths of 0.004-inch-diameter fuel particles ejected in plane 


of rotor, with B = -k5° and 452. Rotor radius (to nozzle), 25 feet; 
tip speed, 1000 fps; internal air speed relative to nozzle, 200 fps; 
fuel-injection velocity relative to nozzle, 250 fps; centrifugal 
acceleration at nozzle, 12408. 
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Figure 11.- Paths of 0, 
of rotor, with! B = 
15 feet; tip Speed, 


Ugal acceleration at 


Radial distance 
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452. Rotor radius (to nozzle), 


100 fps; internal air Speed relative to nozzle, 


fps; fuel-injection velocity relative to nozzle, 150 fps; centrif. 
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Figure 12.~ Paths of 0.002-, 
ejected in plane of rotar, 


D = 0.002 1m.10.004 ind 0,006 in, 


15 


155 


0.00%... and 0.006-inch-diameter fuel particles 
with B = 0°. Rotor radius (to nozzle), 


15 feet; tip speed, 1000 fps, internal air speed relative to nozzle, 
400 fps; fuel-injection velocity relative to nozzle, 150 fps; centrif- 
úgel acceleration at nozzle, 2070g. 


SOVEGI WI VOVN 


GETLO TIS 


LE 


TELOTUISHY 


Angular displacement, y, deg 


245 25 | LAS. ar. 
; Radial diatence fram rotor axie r, ft ] 


Figure 13.- Paths of 0.004-inch-diameter fuel particles ejected in plane 
of rotor, with B = -45° and 45°. Rotor radius (to nozzle), 25 feet; 
tip speed, 800 fps; internal air speed relative to nozzle, 400 fps; 
fuel-injection velocity relative to nozzle, 150 fps; centrifugal 
acceleration at nozzle, 795g. 
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Figure 14... Paths of 0.004-inch-diameter fuel particles ejected in plane 
of rotor, with B = -45° and 45°. Rotor radius (to nozzle}, 25 feet; 
tip speed, 1000 fps; internal air speed relative to nozzle, 400 fps; 
fuel-injection velocity relative to nozzle, 150 fps; centrifugal 
acceleration at nozzle, 12408. 
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Figure 15.- Paths of 0.004-inch-diameter fuel particles ejected in plane 
of rotor, with. 8 = -45°, 0°, and h5?. Rotor radius (to nozzle), 
15 feet; tip speed, 1000 fps; internal air speed relative to nozzle; 
100 fps; fuel-injection velocity relative to nozzle, 250 fps; centrif- 
ugal acceleration at nozzle, 2070g. 
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Figure 16.- Paths of 0. 002-, 0.00h-, and 0.006-inch-diameter fuel particles 
ejected in plane of rotor, with B = 09. Rotor radius (to nozzle), 
lo feet; tip Speed, 1000 fps; internal air speed relative to nozzle, 
400 fps, fuel-injection velocity relative to nozzle, 250 fps; centrif- 
ugal acceleration at nozzle, 2070g. 
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Figure 17.- Paths of O.O0l.inch-diameter fuel particles ejected in plane 
of rotor, with p = -h5? and 45°. Rotor radius (to nozzle), 25 feet; 
tip speed, 800 fps; internal air speed relative to nozzle, 400 fps; 
fuel-injection velocity relative to nozzle, 250 fps; centrifugal 
acceleration at nozzle, 795g. 
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Figure 18.- Paths of 0.004-inch-dismeter fuel particles ejected in plane 
of rotor, with f = -45° and 45%. Rotor radius (to nozzle), 25 feet; 
tip speed, 1000 fps; internal air speed relative to nozzle, 400 fps; 
fuel-injection velocity relative to nozzle, 250 fps; centrifugal 
acceleration at nozzle, 120g. 
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Figure 19.- Path of 0.004-inch-diameter fuel particle ejected at an 


angle of 45? out of the plane of the rotor, with B - 0. Rotor 
radius (to nozzle), 25 feet; tip speed, 800 fps; internal air speed 
relative to nozzle, 200 fps; fuel-injection Velocity relative to 
nozzle, 150 fps; centrifugal acceleration at nozzle, 195g. 
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